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Abstract 

Context. The frequency of planets in binaries is an important issue in the field of extrasolar planet studies, because of its relevance 
in estimating of the global planet population of our Galaxy and the clues it can give to our understanding of planet formation and 
evolution. However, only preliminary estimates are available in the literature. 

Aims. We analyze and compare the frequency of planets in multiple systems to the frequency of planets orbiting single stars. We 
also try to highlight possible connections between the frequency of planets and the orbital parameters of the binaries (such as the 
periastron and mass ratio.) 

Methods. A literature search was performed for binaries and multiple systems among the stars of the sample with uniform planet 
detectability defined by Fischer & Valenti (2005), and 202 of the 850 stars of the sample turned out to be binaries, allowing a statistical 
comparison of the frequency of planets in binaries and single stars and a study of the run of the planet frequency as a function of the 
binary separation. 

Results. We found that the global frequency of planets in the binaries of the sample is not statistically different from that of planets 
in single stars. Even conservatively taking the probable incompleteness of binary detection in our sample into account, we estimate 
that the frequency of planets in binaries can be no more than a factor of three lower than that of planets in single stars. There is no 
significant dependence of planet frequency on the binary separation, except for a lower value of frequency for close binaries. However, 
this is probably not as low as required to explain the presence of planets in close binaries only as the result of modifications of the 
binary orbit after the planet formation. 

Key words. (Stars:) Planetary systems (Stars:) binaries:general (Stars:) binaries:visual (Stars:) binaries: spectroscopic (Stars:) statis- 
tics 



1. Introduction 

The increasing number of extra-solar planets discovered in 
binary or multiple stellar systems (Eggenberger et al. 2004; 
Raghavan et al. 2006; Desidera & Barbieri 2007) suggests that 
planets can form and survive in a variety of stellar environments. 
Determining the frequency of planets in binaries is an important 
issue in the field of extrasolar planets studies, because of its rel- 
evance in estimating the global planet population of our Galaxy 
(more than half of the solar type stars are in binary or multiple 
systems as reported in Duquennoy & Mayor 1991) and because 
of the clues it can give to our understanding of planet formation 
and evolution. The study of the properties of planets in binaries, 
as well as any difference to those of the planets orbiting single 
stars, could shed light on the effects from the presence of the 
companions. 

A recent study by Desidera & Barbieri (2007) shows that 
the mass distribution of short period planets in relatively tight 
binaries (separation < 150 - 200 AU) is significantly different 
from that of planets orbiting the components of wide binaries 
and single stars. There are also other possible peculiar features 
of planets in tight binaries compared to planets orbiting single 
stars, such as a lack of long-period planets and multiple plan- 
ets, that need confirmation. The properties of exoplanets orbit- 
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ing the components of wide binaries are instead compatible with 
those of planets orbiting single stars, except for a possible greater 
abundance of high-eccentricity planets. 

This result implies that the formation and/or migration 
and/or dynamical evolution processes acting in the presence of 
a sufficiently close external perturber are modified with respect 
to single stars. Several scenarios can be devised to explain it, 
such as a different formation mechanism for planets in tight bi- 
naries (e.g. disk instability induced by dynamical perturbations, 
as proposed by Boss 2006), and enhanced migration and ac- 
cumulation rate in the presence of a stellar companions (Kley 
2000), and dynamical interactions after planet formation (Pfahl 
& Muterspaugh 2006). 

However, to better understand the cause of these anomalies 
and the origin itself of the planets in very close binaries (a chal- 
lenge for current planet formation models) a key piece of in- 
formation is missing: the frequency of planets in binaries as a 
function of the binary separation and compared to that of plan- 
ets orbiting single stars. Determining the planet frequency in bi- 
naries is made difficult by the biases against binaries in most of 
the ongoing planet search surveys and by the incompleteness of 
binary and planet detections in these samples. 

A first step in this direction has been made by Eggenberger 
et al. (2006), performing an adaptive-optics search for compan- 
ions around stars with and without planets and without previ- 
ously known stellar companions from the Coralie survey. This 
guarantees a fairly homogeneous binary detectability in their 
sample. However, the small number of objects (and the lack of 



confirmation of the physical association in a few cases at the 
time of the presentation of their preliminary results) did not al- 
low them to make clear inferences on the planet frequency and 
in particular on possible differences as a function of the binary 
separations. 

A much wider sample (850 stars vs the 110 stars studied 
by Eggenberger et al. 2006) that might be used for a study of 
planet frequency is the 'Uniform Detectability' (hereafter UD) 
sample collected by Fischer & Valenti (2005) (hereafter FV05). 
This sample is complete for detecting planets with radial veloc- 
ity (hereafter RV) semi-amplitude > 30 m/s and period < 4 yr. 
However, the binarity of stars in the UD sample has not been 
considered up to now. 

Despite some incompleteness and biases concerning bina- 
rity, this sample can be considered valid to draw an independent 
measurement of the frequency of planets in binary stars, thanks 
to the completeness of planet detection and the large sample size. 
Thus, the binarity of the stars with uniform detectability was in- 
vestigated in this work by searching some stellar catalogs listing 
stellar companions (§ 3). The result is a sub-sample of UD bi- 
naries, separated according to their different values of periastron 
and critical semimajor axis for dynamical stability of planetary 
orbits (see Holman & Wiegert 1999) (§ 3.2). In this way it has 
been possible to compare the values of the frequency of plan- 
ets in the two sub-samples (single stars and binary stars) and to 
verify a possible dependence of the frequency on critical semi- 
major axis and periastron (§ 4). The biases against binaries in 
the original selection sample, the completeness of the binary de- 
tection for stars with and without planets and their impact on the 
results are discussed in § 4.1. The results are discussed further 
in § 5, and in § 6 we summarize our conclusions and offer future 
perspectives. 

2. The uniform detectability sample 

The uniform Detectability sample has been built by considering 
that, despite the detectability of the planets' changes from star 
to the next and from a survey to the other because of the dif- 
ferent time span of the observations and different levels of RV 
errors, we can consider it complete for companions with veloc- 
ity amplitudes K>30 m/s and orbital periods shorter than 4 years. 
Then, beginning from the initial target list, which included 1330 
stars observed by the Lick, Keck, and Anglo Australian Surveys, 
FV05 selected a subsample of 850 stars that satisfy these entries 
provided that at least 10 observations spanning four years were 
available. Stars that were added after a planet was discovered 
by other groups were not included in the sample. However, stars 
independently present in one of these surveys were considered 
even if a planet was detected first by another group. Only plan- 
ets with K> 30 m/s and orbital periods shorter than 4-years were 
considered for the study of planet frequency. This corresponds 
to Saturn-mass planets for the shortest periods and Jupiter-mass 
planets for 4 year orbits. 

2.1. Changes in th UD sample 

During the analysis made for our work, we made some changes 
to the original UD sample, such as: 

- we excluded 5 1 Peg because it was added to the considered 
target lists after the planet detection by Mayor & Queloz 
(1995) (Marcy et al. 1997, Fischer 2005 private communi- 
cations) 



- t Boo and v And were marked in FV05 as "without planets", 
but the known companions fulfill all the selection criteria for 
the UD sample (Fischer 2005, private communication), so 
we included these stars as "with planets"; 

- HD 20782 hosts a planet detected after the compilation of the 
UD sample (Jones et al. 2006), but it is coherent with the UD 
requirements, so that this star has been considered as "with 
planets". The lack of this planet in the original UD sample 
confirms the hypothesis that high eccentricity (e = 0.93 in 
this case) acts to make detection more difficult, as suggested 
by Cumming (2004); 

- HD 18445 is flagged as "with planets" in the UD sample, 
probably because of a typo. It is not listed in the tables of 
stars with planets. The RV companion (Msin; = AAMj) 
has a mass outside the giant planet range, and it was shown 
to be a 0.18 M© star orbiting nearly pole-on by both as- 
trometry (Halbwachs et al. 2000) and direct imaging (Beuzit 
et al. 2004), so we considered it as "without planets" (see 
Appendix A for details); 

- we considered HD 196885 as "without planets" because the 
companion reported in FV05 was not confirmed by Butler 
et al. (2006) (see appendix A for details); 

The modified UD sample that is the result of those changes 
was then searched for companions, in order to build a sample of 
UD binaries. 

3. Searching for UD binaries 

In order to identify known or claimed companions for the stars 
included in the UD sample, we checked available sources list- 
ing stellar companions. Some of the most important sources are 
listed below: 

- The Hipparcos and Thyco Catalogues (Perryman & ESA 
1997); 

- The Catalogue of the Components of Double and Multiple 
Stars (CCDM) (Dommanget & Nys 2002); 

- The Washington Visual Double Star Catalog (WDS)(Wor\ey 
& Douglass 1997); 

- Sixth Catalog of Orbits of Visual Binary Sfars(Hartkopf & 
Mason 2006); 

- The Catalogue of Nearby Stars, Preliminary 3rd Version 
(C3) (Gliese & JahreiG 1991); 

- Gould & Chaname (2004): New HIPPARCOS-based paral- 
laxes for 424 faint stars; 

- Nidever et al. (2002): Radial Velocities for 889 late-type 
stars; 

- Allen et al. (2000): Wide binaries among high-velocity and 
metal-poor stars; 

- Makarov & Kaplan (2005): Proper motion derivatives of bi- 
naries; 

- Tokovinin (1997): MSC - a catalogue of physical multiple 
stars 1 

- Latham et al. (2002) Orbits of 1 71 single-lined spectroscopic 
binaries; 

- Halbwachs (1986) Common proper motions stars in AGK3; 

- The revised NLTT catalogue (Salim & Gould 2003); 

- Valenti & Fischer (2005) Spectroscopic properties of cool 
stars. I. (VF05); 

- The Tycho Double Star Catalogue (TDSC) (Fabricius et al. 
2002); 

1 Updated version available at www.ctio.noao.edu/ atokovin/stars/ 



- A Catalog of Northern Stars with Annual Proper Motions 
Larger than 0.15" (LSPM-NORTH Catalog)" (Lepine & 
Bongiorno 2007) 

- SB9: The ninth catalogue of spectroscopic binary orbits 
(Pourbaix et al. 2004) 

- 2MASS All Sky Catalog of point sources (Cutri et al. 2003) 
has been used only for deriving JHK photometry (used for 
mass determination) and for common proper motion confir- 
mation, not for a search for further companions. 

We also consider additional references for individual objects (see 
Table 8 and Appendix A). 

3.1. Selection criteria 

After this search, we excluded from our UD binary sample 
the stars with non confirmed companions and with companions 
listed in CCDM but with inconsistent proper motions (and with- 
out other indication of binarity found in literature), and consid- 
ered those stars as singles. Stars with long-term RV and/or as- 
trometric trends were included in the binary sample, on the basis 
of the dynamical evidence of a companion. The RV trends we 
included (from Nidever et al. 2002) cause an overall RMS of RV 
faster than 100 m/s that cannot be due to planetary companions. 
We also included stars with brown dwarf companions in the sam- 
ple of binaries. At small separation, the existence of the brown 
dwarf desert (see Butler et al. 2006) guarantees little ambiguity 
in the classification of an object as a massive planet or a brown 
dwarf (a couple of individual cases are listed in Appendix A). 
At large separations, where brown dwarfs companions are prob- 
ably more frequent (Gizis et al. 2001), this ambiguity should 
be taken into account but is again limited to a few individual 
cases, e.g. HD 206860, which has a T dwarf companion of mass 
0.021 + 0.09 M according to Luhman et al. (2007), smaller than 
the mass limit for planetary companions adopted by Butler et al. 
(2006). The small number of brown dwarf companions makes 
this issue irrelevant for the global results. 

3.2. The sub-sample of UD binaries 

The properties of the UD binaries, selected from the modified 
UD sample, are listed in Table 8. The stars with both components 
included in the UD sample are listed twice, otherwise only the 
star under planet scrutiny is listed. If more than one companion 
is known, we report only the closer one, because of its stronger 
influence on planetary formation/stability. For hierarchical triple 
(or higher-order multiplicity) systems, for which the isolated star 
is included in the UD sample, we sum the mass of the closest 
pair to consider its effective dynamical influence. The minimum 
mass is listed for single-lined spectroscopic binaries. 
For each object we report: 

- the HD number; 

- the projected separation; 

- the eccentricity, if the binary orbit is available. Otherwise 
we assume e = 0.31, which is a median value between 
those adopted by Fischer & Marcy (1992) and Duquennoy 
& Mayor (1991). 

- the semi major axis in AU. For those stars for which the bi- 
nary orbit is not available, the semimajor axis was estimated 
from the projected separation a(AU) = 1.31 * p (arcsec) *d 
(pc) (Fischer & Marcy 1992; Duquennoy & Mayor 1991); 
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Figure 1. Values of a crit (solid line) and a cmss (dashed line) and a tii 
(dotted line), versus binary semimajor axis, for e=0.31, M com = 0.5 M Q , 
and M obj =lM Q . 

- the mass of the objects, from Fischer & Valenti (2005) (here- 
after VF05) 2 ; 

- the companion mass 

a) from VF05 if both components are included in the UD 
sample; 

b) calculated with the mass-luminosity relations derived by 
Reid & Gizis (1997); Delfosse et al. (2000); 

c) from other literature sources (listed in the Table caption); 

- the critical semi-major axis for dynamical stability of plane- 
tary companions on circular orbits coplanar with the binary 
orbit, calculated using the equation 

a cr it = (0.464 - 0.380yu - 0.63 \e hin + 0.586yue w „) a bin + 

(0.1504„-0.198^.„) flw „ (1) 

by Holman & Wiegert (1999). InEq. 1: fi = M ^lZ ct>m > a bin is 
the semi major axis and ei, in is the eccentricity of the binary 
orbit. 

We chose reference value, because it is a physical 

quantity that represents, better than the semi-major axis or the 
projected separation, the dynamical effects of the presence of the 
companion first on the circumstellar region and then on planet 
formation and stability. This feature, in fact, includes both the 
orbital parameters and the mass ratio, and represents the max- 
imum value of the semimajor axis for stable planetary orbits 
around the planet hosts. The value of a crit is higher than the limit 
of the region in which the encounter velocities of planetesimal 
is small enough to allow the accretion of kilometer-sized plan- 
etesimals (a cross , Thebault et al. 2006). The radius of tidal trun- 
cation of the circumstellar disk a t a (Pichardo et al. 2005; Pfahl 
& Muterspaugh 2006) is intermediate between a crit and a cross . 
Figure 1 shows the values of a cri , (solid line), a cross , and and, 
versus binary semimajor axis, for fixed values of eccentricity and 
masses (e = 0.3, M com =0.5 M Q , and M ob j-\M Q ). 

4. Results and selection effects 

The result of this search for companions of the UD stars is a sub- 
sample of 202 objects, 15 of those having planets, so the global 

2 the exceptions are listed in Appendix A 



frequency of planets in the UD binary sample is 7.4%. The fre- 
quency of planets in the UD single stars sub-sample is 5.3% (see 
Table 1). The two frequencies are compatible within their er- 
rors 3 , and the slightly higher value of the global frequency in 
the binary sub-sample is probably due to some incompleteness 
in the sample, which is discussed in Sect. 4. 1 . 

The rather large sample size allows us to divide stars in some 
sub-samples according to different values of critical semimajor 
axis for dynamical stability of planets (hereafter a crit ; see Eq. 1). 
All the stars with RV and/or astrometric trend and without direct 
imaging identification or full orbit characterization (37 objects) 
were included in the closest bin, as it is likely that the companion 
responsible for the trend has a small separation. 

We also binned the UD binary sample according with the 
periastron because this allowed us to make a direct compar- 
ison with theoretical expectations, such as those of Pfahl & 
Muterspaugh (2006). The resulting values of the frequency are 
listed in Tables 1 and 2, together with the characteristics of each 
sample. In Table 1 we also show the values of frequency for the 
complete UD binary sample and for the UD single sub-sample. 
Figures 2-3 show a erir and periastron vs mass ratios for the bina- 
ries of the sample with and without planets. 

These figures suggest that the binary components hosting 
planets usually have low-mass secondaries, for both low and 
high separation values. This item is confirmed by Fig. 4 which 
shows the histogram of the mass-ratio values for the two popula- 
tions. In order to better investigate the dependence of the planet 
frequency on the binary mass-ratio, we performed a Kolmogorov 
Smirnov Test (hereafter KST) on the two populations. The re- 
sulting value for the KST probability is ~ 11% for the entire 
sample (excluding objects with trends, for which the value of 
the secondary mass is unknown) and ~ 10% for the stars with 
a cr i t > 20 AU. The difference between the distributions is not 
significant enough to allow us to confirm or reject the hypothe- 
sis that the presence of planets is favored in binaries with a low 
value of the mass ratio. This feature, as discussed in Sect. 4.1.2, 
could also be due to how the stars hosting planets are preferen- 
tially searched for low-mass companions, and this is probably 
one of the causes of their lower observed values of the mass ra- 
tio. On the other hand, the lack of planet detections up to now 
in the ongoing RV planet search using SARG at TNG (Desidera 
et al. 2006), which targets only binaries with similar components 
and with a typical separation between 100 to 400 AU, suggests 
that the mass ratio might be an important parameter for the oc- 
currence of planets. Therefore, the possible role of the mass ratio 
on the frequency of planets remains an open point that requires 
further investigations. 

4.1. Completeness and selection effects 

These results suggest that the planet frequency in binaries and 
single stars is similar. However, selection effects concerning bi- 
naries in the original definitions of the Lick, Keck, and AAT 
samples and possible incompleteness of binary detection might 
alter the results; we discuss them here in order to better estimate 
the frequency of planets in binaries. 



3 The errors reported in Tables 1 and 2 are calculated with the equa- 
tion 
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They do not include the additional error due to the incompleteness of 
binary detection. 
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Figure 2. Critical semimajor axis (from Holman & Wiegert 1999) vs 
mass-ratio for the binaries with planets (filled circles) and without plan- 
ets (open circles) in the UD sample. 
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Figure 3. Periastron vs mass-ratio for the binaries with planets (filled 
circles) and without planets (open circles) in the UD sample. 
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Table 1. Frequency of planets in binaries with different values of a crit . 
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Table 2. Frequency of planets in binaries with different values of peri- 
astron. 
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Figure 4. Distribution of mass-ratio values for star with (solid line) and 
without (dashed line) planets. Panel (a) shows the histogram and panel 
(b) the cumulative distribution. 



4.1 .1 . Estimate of the missing binaries 

The first issue to consider is the completeness of the binary cen- 
sus in the sample. We can derive an upper limit to the incom- 
pleteness by assuming that the stars in the UD sample have the 
same binary frequency as in Duquennoy & Mayor (1991) (57%). 
The number of missing binaries in the sample would then be 
~ 296 4 . 

To derive a lower limit to the frequency of planets in binaries, 
we assume that all the "missing" binaries are without planets. We 
then obtain f bin = 3.09 ± 0.94% and f sin = 9.31 ± 0.69%; i.e. the 
frequency of planets in binaries cannot be less than one third of 
that of planets orbiting single stars. We stress that this is a very 
conservative lower limit on the frequency of planets in binaries 
(an upper limit on the incompleteness of the binary detection in 
the sample), because the UD sample has selection biases against 
binarity. In fact the input target lists exclude spectroscopic bina- 



4 For the estimate of the number of unrecognized binaries, we ex- 
cluded (and then consider them as single stars) the stars with brown 
dwarfs companions, because the substellar companions are not included 
in the statistical analysis by Duquennoy & Mayor (1991); hence, the 
number of binaries found in the UD sample decreases to 188. 



ries and binaries with separations less than 2 arcsec known at the 
time of the target selection (see Jones et al. 2002; Marcy et al. 
2005; Wright et al. 2004). Therefore we expect that the absolute 
binary frequency in the UD sample is significantly lower than 
the unbiased one derived by Duquennoy & Mayor (1991); the 
distribution of the orbital parameters and mass ratio of the bi- 
naries in the UD sample should also be different with respect to 
unbiased samples. 

To take this selection bias in account, we did a Monte Carlo 
simulation making an estimate of the fraction of binaries with 
p <2" expected on the basis of the period and mass-ratio distri- 
bution assumed by Duquennoy & Mayor (1991) 5 . In this way 
we found that ~ 45.7% of the binaries predicted by Duquennoy 
& Mayor (1991) have p < 2" for the distance distribution of the 
stars of the UD sample, so we expect that a significant fraction 
of the 'missing' binaries are given by stars excluded at the time 
of the target selection because of their small separation. 

Combining these features, we find that the total number of 
binaries expected in the UD sample with a separation > 2 arcsec 
is 221. Our census yields 138 binaries (120 without considering 
wide companions orbiting close pairs with a separation < 2 arc- 
sec) so there should be 83 missing binaries (95 without consid- 
ering wide companions in triple systems) with separations > 2 
arcsec. 

These numbers are expected to hold if there are no biases 
besides the exclusion of binaries with separations smaller than 
2 arcsec in the UD sample. However, it is possible that the stars 
that were added later to the samples for specific reasons (e.g. 
high metallicity, Butler et al. 2000; Tinney et al. 2003) have dif- 
ferent selection criteria concerning binarity. In some individual 
cases, the inclusion in the input target lists of certain types of 
systems is favored: a few binary systems with similar compo- 
nents were probably included in the sample after dedicated stud- 
ies of chemical abundances differences between the components 
(Gratton et al. 2001; Martin et al. 2002). Unfortunately, we do 
not have enough information about the details of the building of 
the sample. We do expect that this only plays a minor role in the 
global binary statistics. 



4.1 .2. Completeness of the binarity of planet hosts vs non 
planet hosts 

The role of completeness in binary detection would be minor for 
our purposes if the completeness itself were independent of the 
occurrence of planets. Unfortunately, this is not the case. Planet 
hosts are systematically searched for companions after planet 
discoveries. Therefore, the completeness of the binarity of planet 
hosts is certainly greater than that of stars without planets. This 
bias could cause a spurious increase in planet frequency in bina- 
ries and should be carefully addressed. We proceeded as follows. 
We considered the binarity of planet hosts in the UD sample (in- 



5 The distribution in orbital period assumed by Duquennoy & Mayor 
(1991) is 
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where logP = 4.8, <x logP = 2.3 and P is in days. For the mass-ratio 
q = M 2 /Mi they assume 
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where /i = 2.3, a q = 0.42 and k = 18 for their sample of G-Dwarfs. 
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Table 3. Frequency of planets in binaries with different values of a cril , 
without considering as a binary the planet-host whose companions were 
discovered thanks to dedicated follow up after planet detection. 



eluding those with K < 30 m/s and/or P > 4 yr; 21 objects over- 
all). If the binarity was found on the basis of dedicated studies af- 
ter planet discoveries, the star was classified as single; if instead 
the binarity was known independently of the planet discovery 
(e.g. inclusion in WDS, CCDM, Hipparcos, detection of astro- 
metric and RV trends, etc.), the star was kept as a binary. This 
results in a change in binary status for 7 stars (4 with UD planets 
and 3 with non-UD planets). In this way we obtained a binary 
sample that is less complete than the original one but without bi- 
ases favoring the binarity of planet hosts. The revised frequency 
of planets in binaries then becomes 11/195=0.056 and the cor- 
responding for single stars 38/654=0.058 (the results for all the 
a cr it bins are reported in Table 3). The difference with respect to 
the full sample indicates that indeed the dedicated searches for 
companions around planet hosts somewhat alters the results but 
the frequencies of planets in binaries and in single stars remain 
similar. 

4.1.3. Dependence on separation 

The completeness of binarity in this sample is probably a func- 
tion of the separation. At small separations, the inclusion of stars 
with RV and an astrometric trend probably allows a fairly high 
completeness level. In fact, most of the binaries recently discov- 
ered by means of deep adaptive optics imaging (e.g. HD 13445; 
HD 161797, HD 190406, HD 196885; see App. A for refer- 
ences) would have been included as binaries in this study thanks 
to their dynamical signatures, even without the direct imaging 
identification. Pairs with a small magnitude difference in sepa- 
ration between 0.2 to 10 arcsec should have been detected by 
Hipparcos (Quist & Lindegren 2000). Wide binaries (p > 5 - 10 
arcsec) are more easily discovered and then included in CCDM 
and WDS even for larger magnitude differences. Intermediate 
values of separation (e.g. 1 to 5 arcsec) are probably the most 
incomplete ones, as the detection of a low-mass companion re- 
quires dedicated high-resolution imaging that is not available for 
all the stars, and the companions are not expected to produce de- 
tectable RV or astrometric signatures. 

4.1.4. Effects of incomplete information on orbit and masses 
of the companions 

The determination of a cr i t requires the availability of the full bi- 
nary orbit and an estimate of the mass the companion. All these 
quantities are unknown for the 37 stars with astrometric and/or 
RV trends . In any case, we included these objects in the binary 
statistic because of the significance high level of the trends as 
resulting from the original works by Makarov & Kaplan (2005) 
and Nidever et al. (2002). Furthermore, several objects present 
both the astrometric and RV signatures, which further confirms 
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Table 4. Frequency of planets in binaries with different values of a, 
for the volume-limited sample. 
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Table 5. Frequency of planets in binaries with different values of peri- 
astron for the volume limited sample. 



of their nature as multiple objects. From the timescales of the 
detected orbital motion, we can reasonably infer a period short 
enough to allow the inclusion of these stars in our closest bin 
(dent < 20 AU). However, any more detail on the separation 
distribution of these companions cannot be determined, which 
representing a major limitation for studying of the details run of 
planet frequency vs a crit at small separations. 

For the stars where only the projected separation is avail- 
able, we have considered an approximate values for the semima- 
jor axis and eccentricity values from the adopted statistical rela- 
tions. Even if for individual objects this approximation could be 
quite different from the real value, we expect that statistically it 
would give a realistic representation of the true semimajor axis 
and periastron distributions. 

4.2. The volume-limited sample 

As a further investigation of the role of the biases mentioned in 
Sect. 4. 1 in our results, we have selected a volume-limited sam- 
ple (hereafter VLUD), including the stars within a distance of 
18 pc from the Sun, analogous with FV05. The selection effects 
considered in Sect. 4. 1 are expected to be smaller for the VLUD 
sample. In fact, the census of companions around stars in this 
sample is expected to be more complete than that of the global 
sample because the closest stars were in general more carefully 
searched for stellar companions, and at close distance the 2 arc- 
sec limit corresponds to a smaller physical separation (36 AU 
at 18 pc). The small distance also favors overlap between dif- 
ferent detection techniques. Indeed, only 1 star with astrometric 
trend and without direct imaging identification is included in the 
VLUD sample. 

The volume-limited UD sample includes 129 stars, of which 
44 are known binaries. The fraction of stars with planets in the 
VLUD sub-sample is 9. 1 % for binaries and 9.4% for single stars. 
We selected some bins in a cr „ and periastron also for the VLUD 
sample but, because of the small number of stars included in this 
sample, we considered only 2 ranges of values. Tables 4 and 5 
show the frequency values obtained for the VLUD sub-samples. 
Statistical error bars are larger for the VLUD sample because 
of the small number of objects but the frequency of planets in 
single stars and in binaries is again similar. 

A possibly interesting difference with respect to the full 
UD sample is that the frequency of planets in tight binaries 
(a cr i t < 20 AU) is not lower than in wide binaries and single 
stars. This might be explained looking at Fig. 5. The panels 
show the distribution of a crit for the lowest bin in a cri , for the 
complete UD sample (upper panel) and for the VLUD sample 
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Figure 5. Distribution of a cri , < 20 AU for the complete UD binary sam- 
ple (a) and for the VLUD sample (b). In both panels, the light grey 
column contains the stars with-long term RV or astrometric trends, ar- 
bitrarily placed at a crit < for display purposes. The dark grey column 
represents the contribute of spectroscopic binaries. 



(lower panel). In both panels, the first column on the left con- 
tains the stars included as binaries only because of the long-term 
RV or astrometric trends (i.e. mass and separation of the com- 
panions not known). From these histograms we can easily see 
that the percentage of the stars without definite orbital character- 
istics (thus included only on the basis of dynamical signatures) is 
much larger in the complete UD sample (37/89, ~ 42 %) with re- 
spect to the VLUD sample in which just one object, HD 120780, 
is included only thanks to its astrometric trend. At the same time, 
for the systems for which a cnf can be derived, the distribution of 
the a cr i t values for the VLUD sample is centered on the second 
bin (2.5 < a cr t t < 5.0 AU), while in the complete sample there is a 
relatively large number of binaries (at least 20, excluding those 
with only RV or astrometric trends) for which a cril is smaller 
than ~ 2.5 AU, the separation limit corresponding to P-4 yr for 
solar-type systems. Therefore a lower frequency for planets in 
these systems is expected, as only a part of the separation range 
considered here can host planets on stable orbits. 



5.1. Global estimate of frequency of planets in binary stars. 

Using a sample made of binaries with uniform planet detectabil- 
ity, we compared the frequency of planets in binaries to that of 
planets orbiting single stars. Looking at the whole binary sam- 
ple we can conclude that the frequency of planets in binaries is 
not statistically different from planets orbiting single stars. With 
very conservative assumptions on the incompleteness of binary 
detection, we find that the frequency of planets in the binaries of 
the sample cannot be lower by more than a factor of three com- 
pared to planets orbiting single stars. When considering samples 
that are less affected by selection biases such as the VLUD sam- 
ple, or excluding only the binaries detected thanks to dedicated 
follow-up after planet discoveries, that cause a spurious increase 
of the binary fraction of planet hosts, we found that the two fre- 
quencies are indeed very similar. This is qualitatively consistent 
with the preliminary results of Eggenberger et al. (2006). 

An important point to consider is that this conclusion ap- 
plies to the kind of binaries that are included in UD sample, i.e. 
with a separation distribution biased against close binaries, and 
cannot be directly applied to samples with unbiased binary dis- 
tribution. Furthermore, it applies to the kind of planets fulfilling 
the requirements for inclusion in the UD sample, i.e. P < 4 yr 
and K > 30 m/s. It is possible that e.g. the effects of binarity on 
planets in wider orbits are stronger. 

5.2. Dependence on the binary separation 

The size of our sample allowed us to divide it into some sub- 
groups according to the value of the critical semimajor axis for 
the dynamical stability. In this way we can argue that there is no 
significant dependence of the frequency on a cr « (and on the pe- 
riastron) except for companions with a cr „ less than 20 AU (that 
corresponds to a separation < 50-100 AU, depending on the mass 
ratio of the components). This result makes stronger the conclu- 
sion by Desidera & Barbieri (2007), who reported that the pres- 
ence of distant companions (separation > 300-500 AU) does not 
significantly affect the process of planet formation, as the mass 
and period distribution of planets in such wide binaries are sim- 
ilar to those of planets orbiting single stars (Desidera & Barbieri 
2007), and the frequency of planets is also similar (this work). 

What happens at small separations is more crucial to under- 
standing the role of companions in planet formation. Desidera 
& Barbieri (2007) found that the properties of planets in close 
binaries, in particular the mass distribution, are different from 
those orbiting single stars and components of wide binaries 6 . We 
showed in this study indication for a lower frequency of planets 
in tight binaries. 

The frequency of planets in close binaries can be used to 
further investigate how these planets formed and the origin of 
their anomalous properties. Indeed, Pfahl & Muterspaugh (2006) 
shows that knowing the value of the frequency of planets in close 
binaries 7 should allow two alternative formation scenarios to be 
distinguished. A low frequency (about 0.1% but with an uncer- 
tainty of about one order of magnitude, so we can consider 1 % as 
a limit-value) would be compatible with dynamical interactions 
that cause the formation of the tight binary after planet forma- 
tion. 



6 Note that the definition of tight binaries used by Desidera & 
Barbieri (2007) («„,, < 75 AU) is different than the one adopted here. 

7 Defined as those binaries with semi-major axis less than 50 AU 
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Table 6. Probability of observing 2 binaries with planets in a sample of 
89 (resp. 21) binaries in the UD (resp. VLUD) sample with a cri , < 20 
AU, for different planet frequencies. 
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Table 7. Probability of observing 2 binaries with planets in a sample of 
81 (resp. 19) binaries in the UD (resp. VLUD) sample with periastron 
< 50 AU, for different planet frequencies. 



We tested the probability of obtaining the observed number 
of close binaries with planets in the UD and VLUD samples for 
different frequencies of planets using the binomial distribution. 
This leads to being able to confidently exclude (99%, see Tables 
6 & 7) the preferred value in Pfahl & Muterspaugh (2006). The 
observed frequency is marginally compatible only with the up- 
per limit on planet frequency by Pfahl & Muterspaugh (2006) 
(/ ~ 1%) for the UD sample, and hardly compatible for the 
VLUD sample. The nominal probabilities derived here should 
be taken with some caution because of the possible incomplete- 
ness in binary detection at small separations (anyway estimated 
to be small, in particular for the VLUD sample, see above) and 
because the separation distribution is different than the unbiased 
samples. Nevertheless, our relatively high frequency of planets 
in close binaries suggests that the dynamical interaction after 
planet formation is not the unique channel for creating this kind 
of systems. Instead, we can infer that giant planets might form 
in binaries that have a small separation at the time of planet for- 
mation, possibly in a different way from planets around single 
stars (and around components of wide binaries). 

The run of planet frequency at small separations might shed 
more light on the formation mechanism(s) of planets in binaries. 
Desidera & Barbieri (2007) noted a possible paucity of planets 
in binaries with a critical semimajor axis for dynamical stabil- 
ity in the range 10 - 30 AU; only one planet was found in this 
range, while there are 5 planets with a cri , less than 10 AU and 
4 planets with 30 < a crit < 50 AU. A bimodal distribution of 
planet frequency, with a secondary maximum at a cnf ~ 3 - 5 
AU, is suggested by these data, and it would explain the differ- 
ent characteristics of planets in tight binaries as the result of a 
different formation mechanism. However, our work is not able 
to confirm or reject the reality of such a feature, as only 1 star 
in the UD sample has 10 < a crl , < 20 AU and 1 1 have 20 < a crlt 
< 30 AU (none of them with planets), making the lack of planet 
detections insignificant. Furthermore, the closest bin (a cn( < 20 
AU or periastron < 50 AU) includes several stars (37 out of 89) 
for which a direct detection of the companion is missing and 
which are included only on the basis of the astrometric and/or 
RV trends. Without a determination of the physical parameters 
of these companions, the study of the run of planet frequency at 
small separations is not possible. 



6. Conclusions 

After a detailed search for binarity for all stars in the UD sam- 
ple collected by Fischer & Valenti (2005), we compared the fre- 
quency of planets in binaries and single stars. It turns out that the 
two frequencies of planets are fairly similar. Even taking possi- 
ble incompleteness in the binary detection into account in a very 
conservative way, the frequency of planets in the binaries of the 
sample cannot be more than a factor of three lower than that of 
planets orbiting single stars. 

For moderately wide binaries, the frequency of planets is in- 
dependent on separation. Considering the similar mass and pe- 
riod distributions of planets orbiting single stars and components 
of wide binaries, we then concluded that a wide companion plays 
a marginal role in the formation and evolution of giant planets. 

On the other hand, we found a lower frequency of planets in 
close binaries than in single stars and components of wide bina- 
ries. However, this is probably not as low as required to explain 
the occurrence of planets in close binaries only as the result of 
modifications in the binary orbit after planet formation. This, to- 
gether with the differences in the properties of planets in tight bi- 
naries (Desidera & Barbieri 2007), suggests that planets do form 
in tight binaries in spite of the unfavorable conditions, possibly 
in a different way than for planets around single stars. However, 
crucial issues still need clarification. In fact, it is not yet clear if 
the run of the planet frequency when moving to smaller separa- 
tions is characterized by a continuous decrease, by a sharp cut 
off at which the differences on planet frequency characteristics 
suddenly onset, or by a bimodal distribution, with a minimum 
between 10 < a cri , < 30 AU, a relative maximum at a cn - ; ~ 3 - 5 
AU, a further decrease to zero at extremely small separations, as 
the zone for dynamical stability of planets vanishes. 

These open points might be clarified by a detailed character- 
ization of the binaries in current samples of RV surveys (com- 
pleteness of binary detection and, when possible, full determina- 
tion of the orbital elements) and by the completion of dedicated 
surveys searching for planets in binaries (Desidera et al. 2006; 
Eggenberger et al. 2006; Konacki 2005). 
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Table 8. Properties of binaries found in the UD sample: projected separation (arcsec), eccentricity and semi major-axis (when available), masses 
of the object and the companion, and critical semimajor axis for dynamical stability of planets (Holman & Wiegert 1999). For systems for which 
only the projected separation was available (empty spaces in eccentricity column) the semi major axis was derived from the projected separation 
using the relation a(AU)=1.31p(arcsec)ii(pc) (see Fischer et al. 2002; Duquennoy & Mayor 1991). The asterisk in the last column marks systems 
discussed individually in Appendix A. The mass flag indicates the source for the companion mass: a: M CBmp from VF06; b: M comp from Reid & 
Gizis (1997); Delfosse et al. (2000); c: M comp from individual papers (see Reference below). 
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Remarks: P: Stars with planets as in FV05; S: Stars without planets as in FV05; SB: Spectroscopic Binaries; RV: Stars with RV linear trends 
(see Nidever et al. 2002); A/j: Stars with discrepant proper motion in Hipparcos and Thyco II (see Makarov & Kaplan 2005); G: Stars with 
accelerating proper motions in Hipparcos. (see Makarov & Kaplan 2005). 



References: HD 3074: Allen et al. (2000); HD 4614: Worley & Heintz (1983); HD 4747: Nidever et al. (2002); HD 7693: Nordstrom et al. 

(2004) ; HD 10360 - HD 10361: Worley & Heintz (1983); HD 11964: Allen et al. (2000); HD 13445: Lagrange et al. (2006); Desidera & Barbieri 
(2007); HD 13507: Perrier et al. (2003); HD 13531: Metchev (2006); HD 13612: Worley (1967a); Duquennoy & Mayor (1991); HD 16141: 
Desidera & Barbieri (2007); Mugrauer et al. (2004); HD 16160: Allen et al. (2000), Golimowski et al. (2000); HD 16895: Worley & Heintz 
(1983); HD 18445: Duquennoy & Mayor (1991); Halbwachs et al. (2000); Zucker & Mazeh (2001); HD 20782: Desidera & Barbieri (2007); 
HD 23439: Allen et al. (2000); HD 27442: Chauvin et al. (2006); Desidera & Barbieri (2007); HD 29836: Griffin & Gunn (1981); HD 30649: 
Nidever et al. (2002); Dommanget & Nys (2002); HD 31412: Nidever et al. (2002); Gliese & JahreiB (1991); HD 35956: Vogt et al. (2002); HD 
38529: Reffert & Quirrenbach (2006); Desidera & Barbieri (2007); HD 39587: Nidever et al. (2002); HD 43587: Vogt et al. (2002); Duquennoy 
& Mayor (1991); Salim & Gould (2003); HD 64468: Vogt et al. (2002); HD 65430: Nidever et al. (2002); HD 65907: Tokovinin (1997); HD 
72760: Metchev (2006); HD 77407: Mugrauer et al. (2004), Metchev (2006); HD 86728: Gizis et al. (2000); Lepine & Bongiorno (2007); HD 
90839: Gliese & JahreiB (1991); Duquennoy & Mayor (1991); HIP 52940: Nidever et al. (2002); HD 92222: Fabricius et al. (2002), this paper; 
HD 97334: Burgasser et al. (2005); HD 101177: Duquennoy & Mayor (1991); Salim & Gould (2003); HD 111398: Lepine & Bongiorno (2007); 
HD 120066: Allen et al. (2000); Gould & Chaname (2004); HD 120136: Raghavan et al. (2006); Desidera & Barbieri (2007); HD 120237: Allen 
et al. (2000); HD 120780: Makarov & Kaplan (2005); HD 122742: Nidever et al. (2002); HD 128620 - HD 128627:Worley & Heintz (1983); 
HD 131156: Worley & Heintz (1983); Duquennoy & Mayor (1991); HD 131511: Nidever et al. (2002); HD 134440/39: Allen et al. (2000); HD 
135101: Desidera et al. (2004); HD 139323: Worley & Heintz (1983); HD 139477: Lepine & Bongiorno (2007); HD 140913: Nidever et al. 
(2002); HD 146362 B: Worley & Heintz (1983); Tokovinin (1997); HD 150554: Metchev (2006); HD 156274: Worley & Heintz (1983); HD 
161797: Worley & Heintz (1983); Nidever et al. (2002); Wittenmyer et al. (2006); Makarov & Kaplan (2005); HD 167215: Makarov & Kaplan 

(2005) ; HD 169822: Vogt et al. (2002); HD 174457: Nidever et al. (2002); HD 178911 B: Tokovinin et al. (2000); HD 184860: Vogt et al. (2002); 
HD 185395: Lepine & Bongiorno (2007); HD 187691: Duquennoy & Mayor (1991); HD 190360: Allen et al. (2000); Desidera & Barbieri 
(2007); HD 190406: Liu et al. (2002); HD 191408: Allen et al. (2000); HD 195019: Allen et al. (2000); Desidera & Barbieri (2007); HD 197076: 
Duquennoy & Mayor (1991); HD 196885: Chauvin et al. (2006); HD 198387: Makarov & Kaplan (2005); HD 206860: Luhman et al. (2006) HD 
208776: Nidever et al. (2002); HD 213519: Lepine & Bongiorno (2007); HD 219542: Desidera et al. (2004); HD 219834: Tokovinin (1997); HD 
221830: Allen et al. (2000). 



Appendix A: Comments on individual objects 

A. 1 . List of included binaries 

- HD 3651: planet outside the UD limits. The companion is a 
cool brown dwarf. 

- HD 7693: Hipparcos lists a companion at 0.9 arcsec (~ 19 
AU). Furthermore, the star is listed in CCDM as the wide 
companion of HD 7788, with a separation of 319 arcsec 
(~ 6731 AU). The values of parallaxes and proper motion 
in right ascension would suggest a physical bounding, and 
the RV values reported by Nordstrom et al. (2004) seems to 
confirm this hypothesis. But the RV values for HD 7788 has 
large errors, and the proper motion in declination reported 
by Hipparcos are in disagreement. Probably this peculiarity 
could be explained by considering that both HD 7693 and 
HD 7788 are close binaries themselves. In fact Hipparcos 
lists a companion at 0.9 arcsec for the first star, and one at 
5 arcsec = 105 AU, for the last one. Because we are inter- 
ested in the effect of the binarity on the planetary forma- 
tion/evolution, we will take only the closest companion to 
HD 7693 into account . 

- HD 1 1964: planet outside the UD limits. 

- HD 13445 (GL 86): the companion was discovered by Els 
et al. (2001) and it was classified as a brown dwarf, but suc- 
cessive work shows that the secondary is a ~ 0.5 M Q white 
dwarf (Mugrauer & Neuhauser 2005; Lagrange et al. 2006). 
Desidera & Barbieri (2007) described a possible evolution of 
the system during the mass-loss phase of the originally more 
massive star. The white dwarf companion is responsible for 
the observed RV and astrometric trends. 

- HD 13507: suspected for some time of harboring a planet, 
but the later measurements obtained with ELODIE inval- 
idated this interpretation and instead revealed a classical 
spectroscopic binary velocity curve, caused by a low-mass 
(m sin i ~ 50Mj) companion (Perrier et al. 2003). Adaptive 
optics imaging did not directly detect the companion. 

- HD 13612: triple system. AB is a CPM pair. A is an SB2 
(Duquennoy & Mayor 1991) and the disagreement between 
spectral and photometric parallaxes of A and B are probably 
due to its nature. An additional component at 2' .9 is optical 
(Worley 1967b); the component included in the UD sample 
is HD 13612 B and, because the mass of this star is not listed 
in VF05, we derive it according to the Reid & Gizis (1997) 
and Delfosse et al. (2000) mass-luminosity calibration. 

- HD 16141: planet outside the UD limits. 

- HD 16160=GL105: triple system composed of an inner pair 
and another wide companion (p = 164.8 arcsec = 1200 AU, 
mass 0.38 M ). The orbit of the close pair was derived by 
Golimowski et al. (2000) (a = 15 AU e = 0.75). The mass 
of the close companion is about 0.09 M . 

- HD 18143: triple system. Component B at 6.5 arcsec = 149 
AU, component C at 43 arcsec = 985 AU. 

- HD 18445: is a member of a quintuple system (component 
C). Components AB (HD 18455) form a visual binary with 
P=147 yr and a=1.55 arcsec = 40 AU (Worley & Heintz 
1983), C is at 27.2 arcsec ~ 700 AU and is a spectro- 
scopic binary with minimum mass in the brown dwarf range 
(msin; = 0.042 M G according to Zucker & Mazeh 2001). 
Halbwachs et al. (2000) demonstrated that the pair is close 
to face-on and derived a mass of 0.176 M Q . The close com- 
panion was also visually resolved by Beuzit et al. (2004) at 
0.1 arcsec = 2.6 AU. A further CPM companion (D) is at a 
projected separation 5.090 arcsec = 130 AU from C. 



- HD 20782: listed as a binary in the CCDM catalog. The com- 
panion is HD 20781. Desidera & Barbieri (2007) suggest a 
physical association and consider HD 20782 and HD 2078 1 
as a very wide CPM pair. 

- HD 23439: triple system. The secondary is a spectroscopic 
binary (period = 49 days, masses 0.74+0.12 M Q ). 

- HD 27442 (e Ret): a companion at -13 arcsec ~ 237 AU 
is included in WDS and was confirmed by Chauvin et al. 
(2006). HD 27442 B is probably a white dwarf with a mass 
of about 0.6 M (see Desidera & Barbieri 2007). 

- HD 29836: triple system. HD 285970 is a wide companion 
of the star included in the UD sample. HD 285970 is also 
a short-period spectroscopic binary (Griffin & Gunn 1981). 
Another companion to HD 29836 at 100 arcsec is listed in 
CCDM, probably optical. 

- HD 30649: CCDM lists a companion at 3.37 arcsec = 101 
AU, but probably this is not the cause of the linear trend 
reported by Nidever et al. (2002). 

- HD 31412: additional common proper motion companion 
(CNS3, Lepine & Bongiorno 2007) at 22 arcsec = 792 AU; 
it should not be responsible for the observed RV trend. 

- HD 35956: triple system. The primary (the star in the UD 
sample) is a spectroscopic and astrometric binary ( Vogt et al. 
2002). An additional companion (M = 0.44 M ) is at a pro- 
jected separation of 99 arcsec = 2860 AU. 

- HD 38529: it hosts two low-mass companions with projected 
masses of 0.78 Mj for the inner companion and 12.70 Mj for 
the outer companion. This star shows an astrometric motion 
(it is marked with G-Flag in the Hipparcos Catalogue) which 
is probably due to the presence of the outer companion (HD 
38529 C) for which Reffert & Quirrenbach (2006) derived a 
mass of 2>l + J^Mj, clearly into the brown dwarf regime. A 
stellar companion at very wide separation is also present. 
The classification of this star as a two-planet host with a 
wide stellar companion or a single-planet host with a brown 
dwarf and another wide companion is ambiguous (Desidera 
& Barbieri 2007). 

- HD 43587: triple system. The primary (in the UD sample) 
is a spectroscopic binary (Vogt et al. 2002). An additional 
component (M ~ 0.30 M e , p = 95 arcsec = 2860 AU 
was shown to be physically associated (Duquennoy & Mayor 
1991). And other 3 faint companions are listed, probably all 
optical (NLTT); 

- HD 53705- HD 53706: another distant companion (K5, mass 
0.69) at 185 arcsec = 3000 AU. 

- HD 65907: triple system. The secondary, at a projected sep- 
aration of 60 arcsec from the primary, is itself a close visual 
binary (projected separation 37 AU). Individual masses 0.63 
and 0.35 M . 

- HD 72760: Metchev (2006) reports a companion of 0.13 M Q 
at a separation of 0.96 arcsec = 21 AU, which is probably 
responsible of the astrometric trend found by Makarov & 
Kaplan (2005). 

- HD 77407: 0.30 M G companion at 1.6 arcsec = 50 AU sep- 
aration imaged by Calar Alto Adaptive optic system ALFA 
and confirmed as physically bound to the primary by a multi- 
epoch, high-resolution spectrum (Mugrauer et al. 2004) that 
shows a long-term radial velocity trend forHD 77407 A. The 
companion was also confirmed by Metchev (2006). 

- HD 86728: the M dwarf companion is overluminous with an 
high activity level. Gizis et al. (2000) suggest it is itself a 
close binary. 



HD 90839: CPM pair. There is a third companion, HD 
89862, but it is not physical (see Gliese & JahreiG 1991; 
Duquennoy & Mayor 1991). 

HD 92222 A: it is not included in the Hipparcos Catalogue. 
For this reason the mass of this star is not listed in VF05. We 
derived a photometric distance assuming that both compo- 
nents are on the main sequence; and by using the isochrones 
by Girardi et al. (2002) we found d = 90 pc and M A = 1.09 
M Q , M B = 1.05 AT©. 

HIP 52940 - HIP 52942: triple system. The primary is a spec- 
troscopic binary discovered by Nidever et al. (2002). Both 
components are in the UD sample and both masses are de- 
rived according with Reid & Gizis (1997) and Delfosse et al. 
(2000) mass-luminosity calibrations, because VF05 do not 
list any mass value for these objects. 

HD 97334: a close pair of brown dwarfs (separation 1.5 AU, 
total mass = 0.05 M©) is at a projected separation of about 
2000 AU from the primary (Burgasser et al. 2005). 
HD 99492 B: planet outside the UD limits. 
HD 101 177: triple system. The secondary is a spectroscopic 
binary (P=23 d, masses = 0.74 and 0.31 M ). Component C 
is optical as the other two companions listed most likely are 
(NLTT). 

HD 1 14729: planet outside the UD limits. 
HD 120136 = t Boo: the orbital solution by Hale (1994) 
is very preliminary. Another L dwarf companion candidate 
at 42 arcsec = 664 AU has been reported by Pinfield et al. 
(2006). The physical association has yet to be confirmed. 
HD 120780: CCDM lists a companion at 6 arcsec = 98 AU 
(M=0.53 M ), but probably this is not responsible for the 
astrometric signature reported by Makarov & Kaplan (2005). 
HD 128620-HD 128621 (a Cen): Triple system: another low 
mass companion (Proxima Cen) at very wide separation 
(10000 AU). 

HD 131156=^ Boo: visual binary (p = 4.9 arcsec = 33 
AU, P = 152 yr) (see Worley & Heintz 1983). The vi- 
sual companion explains the long-term radial velocity trend 
with a significant curvature detected by Wittenmyer et al. 
(2006). There are two additional components listed in WDS: 
C (V=12.6,p = 66.7 arcsec = 447 AU, from CCDM), which 
is optical (Duquennoy & Mayor 1991), and D (V = 9.6, 
p = 49 arcsec = 328.3 AU). 

HD 131977: quadruple system. The secondary (HD 131976) 
is itself a binary with a = 0.9 AU and individual masses 0.57 
and 0.38 M©. The fourth component is the brown dwarf GL 
570 B at a projected separation of 1500 AU (Burgasser et al. 
2000). 

HD 1 37778 : wide visual pair with a similar component. Only 
the secondary (HD 137778) is included in the UD sample. 
The primary (HD 137763) is a spectroscopic binary with ex- 
treme eccentricity (e = 0.975; Duquennoy et al. 1992). 
HD 139323: hierarchical triple system. The companion HD 
139341 is a visual binary with a — 18 AU and individual 
masses of 0.74 and 0.81 M©. 

HD 145958: another companion candidate has been reported 
at 0.2 arcsec = 5 AU by Tokovinin (1997) but it needs confir- 
mation (no RV variations observed). Preliminary binary orbit 
of the wide pair in WDS (grade 4). 

HD 146362: member of a hierarchical multiple system. The 
star included in the UD sample is orbiting the nearly equal- 
mass double-lined spectroscopic binary HD 146361 (period 
1.1 days, individual masses 1.10 and 1.09 M©). Preliminary 
orbit of the wide pair in WDS. Another CMP companion, a 
faint M3V dwarf, is at 633 arcsec = 14000 AU. 



- HD 161797 (p Herculis): astrometric orbit by Heintz (1994) 
(period 65 years). The companion was also identified using 
adaptive-optics imaging by Turner et al. (2001) and Debes 
et al. (2002). Radial velocity monitoring revealed a long 
term trend with significant curvature (Nidever et al. 2002; 
Wittenmyer et al. 2006). The star is also listed as a A/j. bi- 
nary in Makarov & Kaplan (2005). Another companion, /i 2 
Her B, lies at 34 arcsec = 285 AU and is itself a visual binary 
with period 43.2 yr (WDS). 

- HD 167215: CCDM lists a companion with p = 54". 8 
(M=1.05 M ), but is not expected to be responsible of the 
astrometric signature reported by Makarov & Kaplan (2005). 

- HD 169822: triple system. The star included in the UD sam- 
ple is a spectroscopic binary detected during the Keck sur- 
vey. Vogt et al. (2002) derived a combined spectroscopic and 
astrometric solution (M = 0.30 M©, a = 0.84 AU). HD 
169889, at 608.4 arcsec = 16426.8 AU, is a CPM star. The 
revised distance of HD 169822 by Vogt et al. (2002) (32 pc), 
coupled with the common RVs for the two stars (Nordstrom 
et al. 2004) makes the indication of a physical association 
stronger. 

- HD 17891 1: triple system (Tokovinin et al. 2000). 

- HD 184860: triple system. The star included in the UD sam- 
ple is a spectroscopic binary (Vogt et al. 2002). The projected 
mass of the companion is in the brown dwarf range (32 Mj). 
An additional companion is at 5.0 arcsec = 151 AU. 

- HD 187691: three components listed in CCDM/WDS. AC is 
a CPM pair, B is optical, C is physical (Duquennoy & Mayor 
1991). 

- HD 190067: Turner et al. (2001) listed a companion at 2.86 
arcsec = 55.19 AU, confirmed by Chakraborty et al. (2002) 
that reported m ~ 0.08 - 0.10M©. 

- HD 190360: Planet outside the UD limits. 

- HD 190406 (15 Sge): A brown dwarf companion (p = 0.79 
arcsec = 14 AU) was found with high-resolution imaging 
made using adaptive optics at the Gemini-North and Keck 
telescopes (Liu et al. 2002). The primary shows a long-term 
radial velocity trend that confirms that HD 190406 B is phys- 
ical, with a minimum mass of M = 48 My. 

- HD 196885: Chauvin et al. (2006) identified at about 0.7 
arcsec a relatively bright companion candidate: HD 196885 
B, that is likely to be a late K-dwarf, with a mass of 0.6 
M Q , orbiting HD 196885 A at a projected physical distance 
of 25 AU. HD 196885 was also classified as a A/j. binary 
by Makarov & Kaplan (2005). The astrometric signature is 
probably due to HD 196885 B. A planetary companion was 
also claimed from the Lick RV Survey 8 , but it is not con- 
firmed in the recent work by Butler et al. (2006). The com- 
panion BD+104351 B, listed in WDS and CCDM, is proba- 
bly optical. 

- HD 197076: the companion listed as B in WDS and CCDM 
(m v =11.6 and p = 94 arcsec = 1974 AU) is optical, but the 
really bound companion is component C (p = 125 arcsec = 
2625 AU) (Duquennoy & Mayor 1991). 

- HD 198387: CCDM lists a companion with p = 12.2 = 511 
AU, but probably this is not responsible for the astrometric 
signature reported by Makarov & Kaplan (2005). 

- HD 215578: the mass of this star is not listed in VF05, so 
we derive the value listed in Table 8 according with Reid 
& Gizis (1997) and Delfosse et al. (2000) mass-luminosity 
calibration. 
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- HD 219542: the masses of the components of this star were 
taken from Desidera et al. (2004), as there is some confu- 
sion in the identification of the components in the photom- 
etry used by VF05 to derive stellar masses (for one of the 
components, the joint A+B magnitude is used). 

- HD 219834: hierarchical triple system. The primary is a 
spectroscopic binary (period=6.2 yr, individual masses of 
0.90 and 0.04 M Q ). The secondary is the star included in the 
UD sample and the mass was apparently derived in VF05 
from the magnitude of the much brighter primary. We as- 
sumed the mass given in the on-line version of the MSC 
(Tokovinin 1997). 

A.2. Unconfirmed binaries 

- HD 23249 (6 Eri): belongs to the group of nearest stars and is 
classified as a weakly active and X-ray soft source (Huensch 
et al. 1998). Fisher et al. (1983) tried to detect a periodic vari- 
ation in the photometric data and suggest that 6 Eri could 
be classified as an RS CVn star: that is, a F-G binary star 
having a period < 14 days, with chromospheric activity and 
with a period of rotation synchronized with its orbital period 
(Linsky 1984), giving the star high rotational velocity induc- 
ing strong activity. This contrasts with the low level of activ- 
ity detected by Thevenin et al. (2005) and the lack of radial 
velocity variations, making this classification doubtful. 

- HD 52265, HD 154857, HD 179949: first epoch observa- 
tions by Chauvin et al. (2006) revealed companion candi- 
dates, but the physical association with the planet hosts has 
not yet been confirmed; 

- HD 102158: the presence of a common proper motion 
companion at 1175.2 arcsec was proposed by Lepine & 
Bongiorno (2007). The corresponding projected separation 
(61000 AU) is much larger than any other companion in 
Table 8, so it is not considered here. 

- HD 107213: Lepine & Bongiorno (2007) propose that this 
star form a wide (projected separation 546 arcsec) com- 
mon proper motion pair with BD+28 2103 = HIP 60061. 
However, the discrepant RVs for the two stars (Latham et al. 
2002; Nordstrom et al. 2004) argue against a physical asso- 
ciation. 

- HD 1 17176: an L dwarf companion candidate at 848 AU has 
been reported by Pinfield et al. (2006). The physical associ- 
ation has yet to be confirmed. 

- HD 168443: as for HD 38529, there is an ambiguity about 
the physical classification of the companion as a massive 
planet or brown dwarf (see Desidera & Barbieri 2007). 

- HD 217107: the presence of the companion listed in WDS 
is not confirmed by recent adaptive optics searches. See 
Desidera & Barbieri (2007) for discussion and references. 
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